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ABSTRACT 

The 3D observed velocities of the Large and Small Magellanic Clouds (LMC and SMC) provide an 
opportunity to probe the Galactic potential in the outskirt of the Galactic halo. Based on a canonical 
NFW model of the Galactic potential, Besla et al.(2007) reconstructed LMC and SMC's orbits and 
suggested that they are currently on their first perigalacticon passage about the Galaxy. Motivated 
by several recent revisions of the Sun's motion around the Galactic center, we re-examine the LMC's 
orbital history and show that it depends sensitively on the dark-matter's mass distribution beyond its 
present Galactic distance. Wc utilize results of numerical simulations to consider a range of possible 
structural and evolutionary models for the Galactic potentials. We find that within the theoretical and 
observational uncertainties, it is possible for the LMC to have had multiple perigalacticon passages on 
the Hubble time scale, especially if the Galactic circular velocity at the location of the Sun is greater 
than ~ 228km s _1 . Based on these models, a more accurate determination of the LMC's motion may 
be used to determine the dark matter distribution in the outskirt of the Galactic halo. 



1. INTRODUCTION 

The Large and Small Magellanic Clouds (LMC and 
SMC) are the two most prominent satellite galaxies in 
the Local Group. They are accompanied by the Mag- 
ellanic Stream which extends along a great circle over 
100 degrees across the sky. One end of the Stream joins, 
both in position and radial velocity, onto the bridge be- 
tween the LMC and SMC (Meatheringham et al 1988). A 
continuous velocity gradient leads to a large infall speed 
at the other end of the Stream (Mathcwson, Cleary, & 
Murray 1974). 

The kinematic properties of the LMC, SMC, and the 
Magellanic Stream have motivated many investigations 
on their origin. One class of models are based on the 
assumption that the Stream is the tidal debris of the 
LMC and SMC during their past perigalacticon passages 
(Toomre 1972, Mirabel & Turner 1973, Fujimoto & Sofue 
1976, Lynden-Bell & Lin 1977, Murai & Fujimoto 1980, 
Gardiner et al 1994, Putman et al 1999). Numerical sim- 
ulations with idealized Galactic potentials have led to the 
prediction for the LMC's proper motion to be 0.02" per 
century in a direction such that it leads the Magellanic 
Stream (Lin & Lyndcn-Bcll 1982). A subsequent ground- 
based observational confirmation gave a slightly smaller 
value (Jones et al 1994, Lin et al 1995). Nevertheless, the 
tidal scenario has been challenged by the lack of excess 
halo stars in the direction of the Stream (Majewski et al. 
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2003) since in such a model, stars and gas should have 
been equally stripped from the Magellanic Clouds. Some 
of these issues may be resolved with an alternative sce- 
nario that the gas in the Magellanic Stream was removed 
from the LMC by a ram pressure as it ploughed through 
the hot residual gas in the Galactic halo (Moore & Davis 
1994). The density of the halo gas must be sufficiently 
high to remove atomic gas from the LMC. It also needs to 
be sufficiently tenuous to avoid any drag on the Stream 
and reduction of its infall speed below its observed values 
(Mastropietro et al 2005). 

One possible test to distinguish between these scenar- 
ios is to reconstruct the LMC's orbital history with a 
set of accurate observational data. The proper-motion 
measurements by Kallivayalil et al. (2006a, 2006b) made 
with the Hubble Space Telescope (HST) provided a set 
of accurate 3D Galactocentric velocities for both LMC 
and SMC. After correcting for the motion of the Sun 
(they adopted a Galactic circular velocity, at the Sun's 
location, of V c = 220km s _1 ), they obtained a transverse 
velocity of the LMC to be ~ 367km s _1 which is much 
larger than the inferred Galactic circular velocity at its 
present-day location and its radial velocity which has a 
positive, modest value of ~ 89km s~ x . 

The newly obtained 3D velocity data confirm that the 
LMC has just passed its perigalacticon. They are also 
useful for the reconstruction of the LMC's orbital history 
especially at the time of its previous perigalactic passage. 
But, such a determination depends on the prescription 
for the Galactic potential. For example, Lin & Lynden- 
Bell (1982) adopted an idealized model which is almost 
certain over simplified. Although, the most recently mea- 
sured proper motion is in agreement with their predic- 
tion, to within a second decimal place, LMC's orbit needs 
to be re-examined with a more appropriate Galactic po- 
tential. Using these velocities and a ACDM-motivated 
MW model with the virial radius (R V i r = 258kpc) and 
mass {M vir = 10 12 M Q ), Besla et al.(2007) reconstructed 
the orbital history of the Clouds and suggested that they 
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are on their first passage about the MW. 

The first-passage conclusion obtained by Besla et 
al.(2007), if verified, would invalidate the tidal origin 
of the Magellanic Stream. In a follow-up investigation, 
Besla et al. (2010) suggested that gas in the Magellanic 
Stream was torn from the SMC by the tide of LMC be- 
fore the two Clouds impinged to the apogalacticon. To 
match the observed column density along the Stream, 
they adopted a gas to star (and dark matter) mass ratio 
for the SMC to be an order of magnitude larger than 
those they adopted for the LMC. In order to distinguish 
these alternative scenarios, it is desirable to check the 
robustness of their results. 

We begin with a re-analysis of the kinematic data. At 
the Clouds' present-day position in the sky, a large frac- 
tion of their observed line of sight and proper motion 
speeds are due to the Sun's motion around the Galac- 
tic center. The distance of the Sun from the Galactic 
center and the circular velocity of the local standard of 
rest have been revised recently (Reid et al. 2009a) to be 
Rq = 8.4±0.6fcpc and V c = 254 ± 16km s _1 respectively. 
These latest determinations are based on the reference 
frame set by the maser sources near the Galactic center. 
In §2, we show that this new value of V c modifies the 
correction of the solar motion and LMC's 3D speed. It 
also significantly revises the reconstruction of the Magel- 
lanic Clouds' orbital history (see a similar conclusion by 
Shattow & Loeb 2009). 

But, there are other recent analysis which lead to a 
range of values for V c 's. Based on 18 precisely measured 
Galactic masers, Bovy et al (2009) deduced a lower esti- 
mation for the circular velocity at the Sun's location to 
be V c = 236 ± 11 km s _1 . In another APOGEE analysis, 
Bovy et al (2012) concluded V c = 218± 6 km s^ 1 and 
claimed that V c < 235 km s _1 at > 99% confidence level. 
In light of this range of V c , we will discuss, in §4.2, the 
dependence of the LMC's orbital history on V c within an 
error bar. 

A high value (254 km s _1 ) of V c also indicates that 
the rotation curve of the Milky Way Galaxy is similar 
to that of the Andromeda Galaxy, suggesting they may 
have comparably massive dark matter halos. Their com- 
bined masses can be estimated from the timing argu- 
ment (Kahn & Woltjer, 1959) because these two most 
prominent galaxies in the Local Group arc currently ap- 
proaching each other and they are presumably bound to 
each other. Based on a recent measurement of M31's 
proper motion, van der Marel et al (2012) estimate a 
total mass of the Local Group to be 3 ~ 4 x 10 12 M Q . 
If their kinematic properties are very similar to each 
other, the mass of either the Milky Way or M31 would 
be 1.5 ~ 2 x 10 12 M Q . This mass is mostly in the form of 
a dark matter halo. In order to verify this possibility, we 
re-analyze the space motion of the LMC based on this 
new information. Wc briefly recapitulate the equation of 
motion in §2. 

In §3, wc illustrate that the LMC's orbital history also 
depends sensitively on the mass distribution in the outer 
Galactic halo, over and above the uncertainties intro- 
duced by the magnitude of V c . As an illustration, wc 
prescribe the Galactic potential as $ oc r~ x (0 < A < 1) 
and show that if A ~ 0.1, the LMC would be currently on 
its second passage about the Milky Way. It is useful to 



adopt a more realistic Galactic potential and determine 
the history of the LMC's orbit. 

There are several attempts to reconstruct the Galactic 
potential based on the velocity information of satellite 
galaxies (Peebles 2001, 2010). However these determi- 
nations are highly uncertain. We follow the approach of 
Besla et al (2007) by utilizing an idealized potential based 
on the ACDM simulations. These simulations produce 
to well determined profiles for the density distribution in 
the inner regions of galaxies (Navarro et al. 1997, here- 
after NFW). However, our illustrative model indicates 
that LMC's past orbit sensitively depends on the poten- 
tial in the outer regions of the galaxy. Numerical simula- 
tions show that there are considerable uncertainties and 
dispersion in the mass distribution for outer regions of 
galaxies. 

In order to explore the possible range of the LMC's 
orbital history, we select, in the §4, several sample dark- 
matter halos of galaxies from numerical CDM simula- 
tions. These models have slight variations in the slope of 
the dark matter density distribution. We also take into 
consideration that with the revised circular velocity of 
the Sun's motion, the mass of the Galaxy may need to 
be upgraded to be comparable to that of M31. With the 
simulated profile of the Galactic potential, we show that 
it is possible for the LMC to have had several encounters 
with the Galaxy. Since the mass ratio between the LMC 
and SMC is 10 : 1, we neglect SMC's contribution. In §5, 
we summarize our results and discuss their implications. 

2. ORBITAL PARAMETERS 
2.1. LMC's spacial velocity 
In their investigation, Kallivayalil et al. (2006a) used 
V c , S un = 220km s" 1 as the circular velocity of the Sun. 
(This value is similar to that derived by Bovy et al. 
2012). They followed the method outlined in van der 
Marel et al. (2002), and obtain the LMC's total, radial, 
and azimuthal velocities as: Vlmc — 378 ± 18km s" 1 , 
VLMC,rad = 89 ± 4km s" 1 , and v L MC,azi = 367 ± 18km 
s _1 respectively. Here we adopt the same procedure to 
correct for the solar motion, including the Sun's pecu- 
liar motion relative to the local standard of rest deter- 
mined by Dehnen & Binney (1998) as {U Q ,V Q ,W Q ) = 
(10.0 ± 0.4, 5.2 ± 0.6, 7.2 ± 0A)km/s, and apply the re- 
cently revised V CtSUn ~ 250km s _1 (Reid et al. 2009a) 
to compute the LMC's total, radial, and azimuthal ve- 
locities. We adopt this high value for V CySun to illustrate 
that it can significantly modified the deduced orbital his- 
tory of the Magellanic Clouds (also see, Shattow & Loeb 
2009). For this revised value of the solar motion we find 
values of Vlmc ~ 356km s _1 , VLMC,rad ~ 63km s _1 , and 
VLMCMzi ~ 351km s" 1 respectively. 

2.2. LMC's equation of motion 

All of the new kinematic quantities (derived with the 
newly revised V CtSun are one sigma smaller than their 
previously determined values and they are expected to 
imply a more bound LMC's orbit around the Galaxy. In 
order to verify this expectation, we compute the LMC's 
orbital history. 

Under the influence of the Galactic potential and dy- 
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namical friction, LMC's equation of motion is 



dr 



$(r) 



DM 



M 



LMC 



(1) 



We adopt the idealized Chandrasckhar(1943) formula for 
dynamical friction such that 
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where X=Vlmc/V^o', cr is the velocity dispersion of the 
DM halo and p(r) is the density of the halo at r. For 
the present application, we neglect the effect of an active 
halo and mass loss from the LMC (Fellhauer & Lin 2007). 
The mass of the LMC is about 2 x 1O 1O M . In most pre- 
vious applications, the Coulomb logarithm in the above 
formula is assumed to be a constant A = b max /b m i n (Bin- 
ney and Tremaine 1987). Here we use a prescription in- 
troduced by Hashimoto et al.(2003) in which b m ax is re- 
placed by the position of the satellite and b min = A.8kpc. 
After the specification of the Galactic potential, equation 
(1) can be solved numerically with a symplectic leapfrog 
integration scheme (Springel et al. 2001). We verify the 
numerical accuracy with a range of time steps. For typi- 
cal models, several hundred steps per orbit are more than 
sufficient for numerical convergence. 

3. DEPENDENCE ON POTENTIAL MODEL OF 
THE MILKY WAY 

We first show the sensitive dependence of the LMC's 
orbital period on the Galactic potential. For illustration 
purpose, we adopt a spherically symmetric, static model 
with a power law potential profile for the Galactic po- 
tential. 

3.1. A simple power law profile 

At LMC's Galactocentric distance, the local circu- 
lar velocity is much larger than its radial velocity and 
smaller than its transverse velocity, and the radial veloc- 
ity is positive. These kinematic properties suggest that 
the LMC has just passed its perigalacticon. Since the 
LMC is close to its perigalacticon, its orbital history is 
mainly related to the structure of the Galaxy outside the 
perigalacticon distance of LMC. Therefore, we mainly fo- 
cus on the outer part of the Galaxy in order to estimate 
how the orbital period of the LMC depends on the model 
parameters. 

We assume an idealized power-law potential $ oc 
r' x (0 < A < 1). The gravity felt by the LMC at a 
distance r is 



F(r) 



GM(r p )M LMC ( r 



-A-l 



(3) 



where r p is the perigalacticon distance of the LMC in 
Galactocentric Coordinates and M(r p ) is the total mass 
of the Galaxy inside r p . Integrate this force to get the 
potential of the Galaxy beyond r p : 
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Fig. 1. — The dependence of ratio r a /r p on A with different 
values of A. From below to above the 10 lines corresponds to A from 
0.1 to 1.0 respectively. The two vertical dashed lines represent the 
value span of A from 0.229 to 0.265 corresponding to Vt = 378km/s 
to 351km/s for LMC. 



Under the conservation of energy and angular momen- 
tum, 



1 



1 + e 
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where the 
r p /r a = (1 
the LMC at r p . 



orbital eccentricity is defined such that 
e)/(l + e), v t is the tangential velocity of 
Substitute $(r p ) and define: 



A(e,X) = 
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In the above equation, the magnitude of A is determined 
by the LMC's azimuthal velocity and distance at peri- 
galacticon. For a given A, each value of A implies a set 
of values for e, r a , and orbital period (sec fig 1). For a 
given value of A, r a /r p would decrease rapidly with in- 
crease of A if the LMC's eccentricity is high. It would 
also decrease rapidly with decreasing value of A if the 
value of A is close to the minimum value required for a 
bound orbit. These tendencies imply that, in the lim- 
itation of high eccentricity, a small modification in the 
value of c (due to fractional changes in the values of A 
or A) would lead to very different orbital periods. 

The derivation of equation (6) is based on the assumed 
conservation of energy and angular momentum. Had we 
taken into account the effect of dynamical friction, the 
loss of energy and angular momentum in time would im- 
ply a larger apogalacticon distance and a longer orbital 
period in the past, albeit the magnitude of its contribu- 
tion is modest. The uncertainties of the LMC's trans- 
verse velocity v t and mass distribution of the Galaxy are 
contained in the values of A and A. In the previous sec- 
tion, we have already discussed the potential reduction 
in the value of vt due to an upward revision in the value 
of V c . sun . For a given value of A, the corresponding in- 
crease in A (eq. 6) would reduce the values of r a j r p (Fig. 
1) and LMC's Galactic orbital period. 
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Fig. 2. — The dependence of orbital history on the LMC's velocity 
with A = 0.1. The dotted lines are the results without dynamical 
friction. Both velocities would lead to more than one passages 
during Hubble time even with under the influence of a strong dy- 
namical friction effect. 

Adopting the most recently revised circular velocity 
of the local standard of rest and LMC's location (r p ~ 
50kpc) as vt and r p , equation (1) can be solved for differ- 
ent values of A. In order to illustrate this point, we choose 
A = 0.1 and M(r p ) = 3.8 x 10 n M o (model A 1 of Klypin 
et al. 2002) to examine the orbital history. This value 
of M{r p ) is identical to that used by Besla et al (2007). 
With these parameters the corresponding value of A for 
LMC is 0.229 and 0.65 under the tangential velocity as 
Vt = 378fcm/s and Vt = 351fcm/s respectively(Fig. 1). 
We choose a small value for A because if the rotation 
curve remains relatively flat in the outside part of the 
Galaxy, A would be very close to zero. 

In order to determine the magnitude of dynamical fric- 
tion from equation (2), we adopt the density distribu- 
tion of the Galactic dark matter halo from equation (4). 
We choose the largest value (1) for X. This approxima- 
tion slightly over-estimates the contributions from the 
dynamical friction. With this form of the potential, we 
obtain more than one past pcrigalacticon passage by the 
LMC for both the previous (which yields v t = 378km 
s _1 ) and the latest solar circular velocities (which yields 
356km s" 1 ) (Figure 2). The smaller magnitude of the 
newly determined azimuthal speed decreases the LMC's 
orbital period significantly as expected. 

4. NFW MODEL AND SIMULATION DATA 

4.1. Models of Galactic potential 

In order to construct a more realistic Galactic poten- 
tial, we separate the baryonic and dark matter contribu- 
tion to the Galactic potential such that 



M(r)=M b (r)+M dm (r) 



(7) 



We assume that the dark matter density profile is de- 
scribed by a NFW model (Navarro et al. 1997): 



Phalo{r) 
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In the above equations, C is the halo concentration, M vir 
and r v i r are the halo virial mass and radius. The p cr is 
the critical density of the universe and 8th is the over- 
density of a collapsed object according to the spherical 
top- hat model(Gunn & Gott 1972). The r vir is defined 
as the radius inside which the average density equals to 
the virial density (p vir = 8 th p M = 5 th p cr * Q. M ). 

In order to be consistent of the simulation param- 
eters, we take 8th = 200, Hubble constant Hq = 
72kms~ 1 Mpc~ 1 for a flat universe, and Om = 0.216 for 
our cosmological model. With these standard paramc- 
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(13) 

Only two independent parameters, C and M V i r , need to 
be specified to define the values of all other halo quanti- 
ties. 

The reconstruction of the LMC's orbital parameters 
depends less sensitively on the Galactic mass distribu- 
tion at distances much smaller than r p where most of 
the baryonic matter resides. Nevertheless, we include the 
baryonic matter's contribution to the Galactic potential 
so that our model is self consistent, ie it can reproduce 
the observed circular velocity of the local standard of 
rest. The baryonic components include the mass of cen- 
tral black hole, the Galactic nucleus, bulge and an expo- 
nential disk with scale length r& We follow the equation 
outlined in Klypin et al. (2002) such that 

M b (r)=M BH + 0.025M Mir [l - exp{-2Mr 115 )} 
+0.142M Mir [l - (1 + r^expi-r 1 - 5 )] 



+0.833M& 



The disk surface density is: 



exp 



r 



£(r) = E exp I — - 



14) 



(15) 



where Mb, v ir is the total mass of the cooled baryons, as a 
fraction of the virial mass of dark matter halo, M b ^ V i r = 
0.05M„ ir and = 3.5kpc is the scale length of the disk 
corresponding to the location of the Sun as 8.5kpc. e.g. 
if M vir = 1O 12 M , the total mass of the disk is 4.2 x 
10 10 M Q and the bulge mass as 7.1 x 10 9 M Q . 

Since the pericenter of the LMC's orbit is about 50kpc, 
which is far away from the typical size of disk and bulge, 
the exact baryonic mass distribution is not very impor- 
tant for the orbital history, but it is relevant to the solar 
circular velocity The rotation curve of an exponential 
disk is (Freeman 1970): 



f(x) = ln(l + x) 



(10) 



4*GEo(3l)Mv)K (v) - h{y)Ki (y)] (16) 

where y = r/ra and I n ,K n are the Modified Bcsscl Func- 
tions of the first and second kind(Binney & Tremainc 
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Fig. 4. — The rotation curve for M vir = 1.5 X 1O 12 M , R vir = 
232fcpc and C = 13. The circular velocity at 8.5fcpc is V c ,sun ~ 
257kms~ 1 . 

2008). The zero-pressure rotation velocity can be deter- 
mined as: 



V c 2 (r) 



disk 



G 

(r) + -[M BH + M nu i +bul (r) + M halo (r)} 



(17) 

Here we do not consider the adiabatic contraction of the 
dark matter halos. 

With this velocity of the Sun, we can calculate the 
current velocity of LMC. If we choose M V i r to be of the 
order 10 12 M Q or 1.5 x 10 12 M Q , it self-consistently gives 
V c of the Sun to be 220km s _1 or 250km s _1 respectively 
from the rotation curve (see Figs 3 and 4). 

When we choose different parameters of NFW models 
to calculate the orbital history of LMC, the V c will also be 
adjusted slightly according to the changing dark matter 
mass inside the location of the Sun. 

4.2. Simulated models of dark matter halo 

The NFW model is well suited for the part of the 
Galaxy within the virial radius. If we extrapolate this 
model for the entire dark matter halo, its density would 
decrease as pocr -3 for r >> r s . However, from the nu- 
merical simulations carried out by Oser et al.(2010), we 
find that the actual halo density profile is somewhat flat- 
ter than expected for the NFW model, ie there is more 



residual dark matter at large distance than that inferred 
from the NFW model. As we have discussed above, the 
LMC's orbital history is very sensitively determined by 
the density gradient of the outside part of the Galaxy. 
Here, we use a more realistic dark halo profiles which is 
extracted directly from the cosmological simulations. 

In order to constrain the dark matter density profiles 
at large Galactocentric radii, we choose a set of 4 'zoom- 
in' cosmological simulations of individual halos (Oser et 
al. 2010). The initial condition of this simulation is 
based on a flat WMAP3 cosmology (Spergel et al. 2007) 
with model parameters: h=0.72, fij, = 0.044, fi m = 
0.216, fl A = 0.74 and er 8 = 0.77. The initial slope of 
the power spectrum is n s = 0.95. 

At redshift zero we identify halos with the help of a 
friends-of-friends finder and choose relaxed halos with 
no massive companion close-by for re-simulation. We 
trace back in time all particles closer than two times the 
virial radius to the halo center in any given snapshot and 
replace those particles with multiple dark matter parti- 
cles of lower mass while adding the small scale density 
fluctuations with the help of GRAFIC2 (Bertschingcr 
2001) The dark matter particles outside the region of 
interest are merged (depending on their distance to the 
re-simulated halo) to reduce the particle count and the 
simulation time. This way we are able to simulate struc- 
ture formation in the cosmological context at high reso- 
lution in a reasonable amount of time. To evolve the high 
resolution initial conditions from redshift z = 43 to the 
present day we use the parallel Tree-code GADGET-2 
(Springel 2005). 

The set of 4 halos in Table 1 have masses in the range 
of 0.5 — 2 x 10 12 M Q /i _1 . The computational domain cor- 
responds to a cosmic cube with the size of lOOMpc/h. 
The dark matter is represented by particles with a mass 
2.5 x 10 7 Mq/i -1 and a comoving gravitational softening 
length of 890/1" 1 pc. The radial dark halo density struc- 
ture is well resolved in the region from r p = 50 kpc to 
the virial radius and beyond. 

In order to apply the simulation data to the compu- 
tation of the LMC's orbit, we first match the simulated 
dark matter distribution inside 200 kpc with NFW mod- 
els to extract the best fitted values of p s and r s . Based on 
these two values, we then match the density distribution 
at around 200 kpc with the fitting formula 



Ps 



{r/r s )(l + r/r s y 



(18) 



to determine the value of /3. Instead of the idealized 
NFW profile in equation (8) which is equivalent to the 
case of (3 = 2, we use the above fitting procedures to 
approximate the dark matter halo's density profile be- 
yond 200 kpc (see Table 1). In Figure 5, we show 
the LMC's orbital history for various model parameters. 
Since the M V i r values for these halos are relatively high 
and j3 < 2, it is possible for LMC to have more than one 
peri-galacticon passage even for V sun < 250km/s. Note 
that there is one low-Myj r model (1167), where the LMC 
is undergoing its first peri galacticon passage. 

In the results listed in Table 1, we also calculate the 
theoretical circular velocity of the Sun, in accordance 
with the Equations (16) and (17). The actual observed 
circular velocity at the Sun's location is related to the 
precise mass structure of the baryonic mass, because for 
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Fig. 5. — Orbital history of the LMC for simulated dark matter 
halos. The dark matter structure beyond 200fcpc is determined 
from equation (18). The solid lines show the orbits for dark halos 
with /3 as given in Table 1. The dotted curves show the situation 
if one would take a dark halo with similar virial parameters but 
fixing beta at 2 as expected for an idealized NFW profile. 



TABLE 1 



data 


M vir (10 12 M @ ) 


r v ir(kpc) 


r a (kpc) 




V c (km s" 1 ) 


0959 


1.67 


240 


21 


1.78 


247 


1061 


1.43 


228 


18 


1.61 


224 


1071 


1.47 


231 


21 


1.69 


230 


1167 


1.29 


221 


21 


1.83 


220 



From the simulation results of Oser, et al. (2010). 

the inner part of the Galaxy, the contribution of circular 
velocity from the baryonic matter in the disk is compara- 
ble to that from dark matter (Fig 3 & 4). The theoretical 
model wc used here for an exponential disk may lead to 
an underestimation of the Sun's circular velocity. Al- 
though the theoretical circular velocities for the models 
listed in Table 1 do not match exactly the latest observa- 
tional data (ie V c ~ 254/cm/s), the extra mass we need to 
make the velocity to be consistent with the observation 
is very small compared with the total mass within 50kpc. 
Since the LMC does not venture into this inner region of 
the Galaxy during its orbital history, the precise struc- 
ture inside 8.5kpc does not affect our determination of 
the LMC's orbital history. 

The obscrvationally determined values of V c remain un- 
certain within the range between 220 km s" 1 to 250 km 
s" 1 (Reid et al. 2009a, Bovy et al. 2009, 2012, also 
see §1). If these values are applied to the NFW models, 
the lower limit of V c would imply a single pcrigalacti- 
con passage for the LMC whereas with the upper limit 
of V c , the LMC would have had multiple perigalacticon 
passages. However, with the modified density distribu- 
tion for the outer galactic halo, the LMC in three mod- 
els (0959, 1061, and 1071) with V c > 228 km s" 1 have 
had multiple perigalacticon passages. Only in one model 
(1167) with V c = 221 km s -1 , the LMC's orbital period 
is marginally smaller than 13 Gyr. Based on these re- 
sults, we suggest if V c > 228 km s _1 , LMC would likely 
have undergone more than one perigalacticon passages. 

The correction resulting from the most recent observed 




Time (year) x1() a 

Fig. 6. — Orbital history of the LMC with the minimum trans- 
verse velocity Vt = 348 km s — 1 . This value of Vt is within the 
error-bar from that (Vt = 367km s" 1 ) derived with V c = 220km 
s -1 . For both cases, we used the same current radial velocity for 
the LMC as V r = 89km s — 1 and the halo potential in accordance 
to model 1167. 

circular velocity of the Sun leads to a smaller transverse 
velocity, orbital period, and apo galacticon distance for 
the LMC (see §2). However, even without this change in 
the circular velocity of the Sun, the la error bar in the 
LMC's observed velocity is ~ 18km s _1 . In order to as- 
sess the implication of these uncertainties on the LMC's 
orbital history, we use a range of transverse velocities 
of the LMC's current orbit that are consistent with the 
uncertainties in Vt of367 ± 18km s _1 and re-calculate 
the orbital history with the simulated halo based from 
the model 1167 which has M mr ~ 1.3 x 1O 12 M and 
R V i r ~ 220kpc. In Figure 6, we show that the velocity 
uncertainties are large enough to introduce a significant 
modification in the LMC's orbital period as we have al- 
ready suggested in §3.2. Although the mean and high 
values of the transverse velocity imply that the LMC has 
just passed through its perigalacton for the first time, 
the lowest value of the transverse velocity would imply a 
multiple peri galacticon passage for the LMC. This error 
analysis poses a challenge to the statistical significance of 
the single passage result obtained by Besla et al (2007). 

4.3. Model including halo evolution 

Since the LMC's orbital period is either comparable 
(in the case of a single passage) or a significant fraction 
(in the case of multiple passages) of the Hubble time, 
it is relevant to consider the evolution of the Galactic 
dark matter halo. However, the build up of the dark 
matter halo is a very complex process. It involves merger 
events, accretion and dynamical relaxation. In order to 
construct a model for the evolution of the Galactic halo, 
we need to identify the dominant mechanism. 

We begin with an observational interpretation. The 
age of the thin Galactic disk can be traced by the age 
(> 10 Gyr) of old white dwarfs in it(Hansen et al. 2004). 
The preservation of the thin disk structure for this pop- 
ulation of stellar remnants suggests that the Galaxy may 
have not had any major merger events during the past 
10 Gyr. However, a significant mass increase due to a 
smooth infall or cold streams cannot be ruled out. 

Another approach to model the evolution of the dark 
matter halo is to utilize the ACDM models which 
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Fig. 7. — Mass evolution according to eq.(19). 



were used to model the density distribution we used in 
§4.1. Here we adopt a recently developed prescription 
(Krumholz & Dckcl 2012) based on the assumption that 
the in-streaming of baryonic and dark matter has led to 
the growth in the Galactic potential and the virial mass 
at a rate 

Mvir,X2 = -0.628Mi ir 14 i2^ (19) 



where M vir V i = M vir /10 1 M Q . In the above equation, 
uj = l.Q8/D(t) is the self-similar time variable of the 
extended Press-Schechter (EPS) formalism, and D(t) is 
the linear fluctuations' growth function. Based on this 
EPS formalism, Neistein et al. (2006) and Neistein & 
Dekel (2008) estimate 



-0.0476[1 



0.093(1 + z) 



-1.22)2. 



5 Gyr _1 . (20) 



The virial radius is related to both virial mass and 
redshift via r vir ocM^ / H(z)' 2 / 3 (Burkcrt ct al. 2010), 
where H(z) is: 



H(z) = H Q [n A + (l - fi A - fl M )(l + z f + M (1 + zf 



The redshift is also related to universal time through: 



(21) 



(22) 



With these two relations, we can model the time evolu- 
tion of the virial radius and calculate the orbital history 
of the LMC with this simplified evolution model for the 
dark matter halo. 

We consider a low-mass (with the virial mass M V i r 
I x 10 12 M Q ) and a high-mass (M mr = 1.66 x 10 12 M Q ) 
case. These limits are comparable with the minimum and 
maximum masses of the simulation data. For both cases, 
we set the concentration parameter C = 12 and neglect 
its evolution. For these two sets of model parameters, the 
initial mass accretion rate is relatively high ( > 100M© 
yr" 1 ) at redshift z > 2(Fig.7). Although in the evo- 
lutionary models, M v i r increased steadily LMC's orbits 
since the previous perigalacticon passage for the evolu- 
tionary models are similar to those with a fixed halo. 
For the low mass model (Fig. 8), the previous apo and 
perigalacticons occurred at 4 and 8 Gyr ago when M V i r 
has already attained 85% and 50% of its value at z = 



no evolution 
with evolution 




Time (Year) 



Fig. 8. — Orbital history of the LMC with and without evolution 
of the dark matter halo. The concentration parameter is C = 12 as 
constant with the current M V i T = 1.3 X 10 Mq and kept constant 
at a value of R v i r = 22lkpc. The corresponding circular velocity 
is Vsun = 236km s" 1 . 

respectively. The late mass increase of the halo poten- 
tial does not modify LMC's orbit significantly. However 
more than lOGyr ago, M v i r was substantially smaller 
than today. The play back integration of LMC's orbits 
for the evolving and the non evolving halo models there- 
fore diverges at high redshift. This divergence is similar 
to that between models which include or neglect contri- 
butions from dynamical friction. However the dynamical 
friction effect decreased at early epochs in the evolving 
halo model due to a decreased density of dark matter 
back then. 

Similar results are also obtained for the high mass 
model (Fig. 9). In this case, the previous apo and peri- 
galacticons occurred at 1.8 and 3.6 Gyr ago when M V i r 
has already attained 94% and 88% of its value at z = 0. 
Thus, the divergence between the orbits for the evolv- 
ing and non evolving dark matter halo is only significant 
more than 8 Gyr ago, when the LMC passed through two 
peri galacticon passages prior to the present epoch. We 
note that when the evolution of the halo is taken into 
account, the number of passages declines by at most one 
in comparison with that for the non evolving halo. 

In the above models, we adopt a constant concentra- 
tion parameter. For an alternative possibility, we con- 
sider a model in which C increased linearly from C = 4 
to C = 12 during cosmic time for the evolving halo mod- 
els. We consider the case with the asymptotic virial mass 
M vir = 1.66*10 10 M Q (see Fig.10). The small concentra- 
tion at early time makes the LMC's orbit more loosely 
bound. However, in comparison with the non-evolving 
halo model (with a fixed C = 12), the difference results 
from C variations is not very significant. It seems that 
the orbital period of LMC is more dependent on the cur- 
rent concentration of the dark matter halo than the ini- 
tial value. 

5. SUMMARY 

Besla et al (2007) utilized LMC's three dimensional 
velocity data to reconstruct LMC's orbital history. They 
made a bold suggestion that the LMC has passed the 
perigalacticon of its Galactic orbit for the first time. This 
conclusion, if it can be verified, would invalidate the tidal 
disruption hypothesis for the Magellanic Stream. 



Zhang, Lin, Burkert & Oser 
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Fig. 9. — Orbital history of the LMC with and without the evo- 
lution of a larger dark matter halo. The concentration parameter 
is C = 12 and constant with the current M v i r = 1.66 * lO 12 Af0 
and Rvi r = 2A0kpc. The corresponding circular velocity is V BU n = 
259km s" 1 . 
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Fig. 10. — Orbital history of LMC with linear evolution of con- 
centration parameter of DM halo, compared to the results with 
fixed C = 12 and C = 4. The final mass is M vir = 1.66 X 1O 12 M . 

We show that these previous results depend critically 
on the LMC's transverse velocity and the dark matter 
distribution in the outer Galactic halo. At the LMC's 
latitude and longitude in the sky, a large fraction of its 
observed proper motion is due to the Sun's circular veloc- 
ity around the Galactic center. We use the most recently 
measured solar circular velocity (254 km s _1 ) rather than 
its conventional value (220 km s _1 ) to deduce a smaller 
transverse speed, apogalacticon distance, and orbital pe- 
riod for the LMC. It raises the possibility that the LMC 
actually had two or more previous pcrigalacticon pas- 
sages. 

We adopt the current velocity of LMC with a differ- 
ent solar velocity, and a simply modeled potential profile 
of the Milky Way outside 50kpc to examine its orbital 
history. We find that the number of passages depends 
sensitively on the slope of the potential and density dis- 
tribution of the dark matter halo as well as the current 
transverse velocity of LMC near its perigalacticon dis- 
tance. 

With an illustrative model (in Fig. 2), we show that, 
for a relatively fiat Galactic potential (with lambda=0.1), 
a modest (7 percent) decrease in the deduced value of V t 



can lead to a much larger (30 percent) decline in the 
LMC's inferred orbital period. Thus, a precise determi- 
nation of LMC's orbital history can provide a sensitive 
measurement on the distribution of dark matter halo in 
the extended outer regions of the Milky Way Galaxy. 

In this paper, we have neglected the perturbation on 
the LMC from any other member of Local Group, in- 
cluding M31. If we choose a massive model of dark mat- 
ter halo (M V i r = 1.66 x 10 12 M Q ), the pericenter's loca- 
tion of the last two orbital periods would be well within 
400kpc (Fig. 9). The current distance between M31 and 
the Milky Way is about 700kpc, it was even further away 
in the past since they are currently approaching each 
other. The M31-Galaxy orbital plane also appears to be 
inclined to the LMC-Galaxy orbital plane. It's pertur- 
bation on the previous the LMC's orbital history appear 
to be weak. This assumption is also consistence from the 
lack of warp along the great circle which is traced out by 
the Magellanic Stream. Nevertheless, a more comprehen- 
sive set of simulations is needed to justify the dynamic 
independence between M31 and the Milky Way galaxies. 
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